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The results of mid-rapidity (0 < y < 0.8) neutral pion spectra over an extended transverse mo- 
mentum range (1 < pr < 12 GeV/c) in y/s NN — 200 GeV Au+Au collisions, measured by the STAR 
experiment, are presented. The neutral pions are reconstructed from photons measured either by the 
STAR Barrel Electro-Magnetic Calorimeter (BEMC) or by the Time Projection Chamber (TPC) 
via tracking of conversion electron-positron pairs. Our measurements are compared to previously 
published n ± and ty° results. The nuclear modification factors Rct> and Raa of n° are also pre- 
sented as a function of pr ■ In the most central Au+Au collisions, the binary collision scaled tt° 
yield at high pr is suppressed by a factor of about 5 compared to the expectation from the yield of 
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p+p collisions. Such a large suppression is in agreement with previous observations for light quark 
mesons and is consistent with the scenario that partons suffer considerable energy loss in the dense 
medium formed in central nucleus-nucleus collisions at RHIC. 

PACS numbers: 25.75.Dw, 13.85.Ni 



I. INTRODUCTION 

The observation of "jet quenching" [1, 2] in central 
Au+Au collisions is one of the most exciting experi- 
mental discoveries at the Relativistic Heavy-Ion Collider 
(RHIC). Experimental signature of this observation in- 
cludes the suppression of inclusive hadron yields at high 
transverse momentum (px) [3, 4] and of associated pr >2 
GeV/c particles on the away-side of a high pr trigger 
hadron [5] . These measurements indicate that RHIC has 
produced high energy density matter that is opaque to 
high pr quarks and gluons [1]. Theoretical calculations 
based on energy loss of high pr partons through gluon ra- 
diation can explain the suppression of light quark mesons 
[6] . Measurements of 7r° at high px provide a fundamen- 
tal tool for probing the parton energy loss in the medium 
created in central nucleus- nucleus collisions at RHIC. On 
the other hand, this medium appears to be transparent to 
direct photons, of which the nuclear modification factor 
Raa is found to be approximately unity at high pr [7] . 
Measurement of the ir° spectrum over an extended pr 
range is a prerequisite to understand the decay photon 
background of the direct photon analysis. This measure- 
ment also provides an important cross-check for other 
pion measurements at RHIC using different detectors. 

In this article we present the first results for the 7T° 
spectra and nuclear modification factors at mid-rapidity, 
over a broad pr region (1 < pr < 12 GeV/c) in Au+Au 
collisions at ^snn — 200 GeV measured by the STAR 
experiment. Neutral pions are reconstructed via the di- 
photon decay channel. Only the west half of the STAR 
Barrel Electro-Magnetic Calorimeter (BEMC) [8] was 
completed and commissioned to take heavy-ion collision 
data in 2004. Measurement of a 7r° spectrum under the 
high multiplicity environment in central Au+Au colli- 
sions is challenging due to the large transverse size of the 
STAR BEMC towers (0.05x0.05 in Ar/x A0), resulting in 
high occupancy and appreciable background contamina- 
tion. The BEMC provides STAR with a trigger capabil- 
ity on high pt photons based on large energy deposition 
in a single BEMC tower or a tower patch. These trig- 
gered BEMC photons can be used to reliably construct 
7r mesons in the high pt region. However, this is not 
possible in the low px region as the energy resolution of 
the BEMC is not sufficiently good. The STAR Time Pro- 
jection Chamber (TPC) [9] has been used to reconstruct 
photons that convert to electron-positron pairs [10]. Ex- 
cellent detection resolution on 7r° invariant masses has 
been achieved from TPC conversion photons. However, 
the small photon conversion probability in the STAR 
detector system restricts the pr reach. By combining 
BEMC photons from high pr triggers and conversion 



photons from the TPC, we have been able to achieve 
good invariant mass resolution on ir° reconstruction and 
measure its spectrum over a broad range of pr- 

II. DATA ANALYSIS 

A. Data Set 

The data used in this analysis were taken during the 
2004 RHIC run for Au+Au collisions at the energy y / s]v"/v 
= 200 GeV. The primary STAR detectors used for this 
analysis arc the TPC and BEMC. A Barrel Shower 
Maximum Detector (BSMD) [8] at a depth of 5 radia- 
tion lengths (Ao) inside the BEMC measures transverse 
shower shape and position with higher precision than 
the BEMC tower. Three Au+Au data sets were used: 
11 x 10 6 events selected by a Minimum-Bias trigger (MB), 
17 x 10 6 events selected by a central trigger, and 2.4 x 10 6 
events selected by a High Tower trigger (HT) . The central 
trigger corresponds to the highest 12% charged particle 
multiplicity events as determined by the coincidence of 
the Central Trigger Barrel and the Zero Degree Calorime- 
ters [11]. The HT trigger, which depends on pseudo- 
rapidity, requires that at least one BEMC tower has de- 
posited transverse energy greater than the HT energy 
threshold of 3-4 GeV. The HT trigger enhances selec- 
tion of events containing high pt photons, and thus helps 
to extend our measurement to higher pt- More details 
about the STAR trigger system and trigger configuration 
can be found in Rcf. [11]. In this analysis, the position of 
the primary vertex is required to be within ±20 cm of the 
center of the STAR TPC along the beam line. This re- 
quirement restricts our conversion photon candidates to 
mid-rapidity, where the detector geometry is relatively 
simple and the material is well studied for reconstructing 
conversion photons. 



B. Photon Identification 

There are two ways to identify photons in STAR: 
The STAR BEMC and BSMD measure photons directly 
from the electromagnetic shower (EMC photon); or the 
STAR TPC reconstructs photon conversion to e+/e~ 
pairs (TPC photon) in materials such as the beam pipe, 
the Silicon Vertex Tracker (SVT), the Silicon Strip De- 
tector (SSD), and TPC walls and gas. In total these 
materials are estimated to be equivalent to about 0.1A , 
with a 10% uncertainty based on studies of conversion 
probability correction. The photon conversion probabil- 
ity will be discussed in Sec. II D and HE. 
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FIG. 1: The di-photon invariant mass distributions using the EMC-TPC method in 0-20% Au+Au collisions at ^snn = 200 
GeV. The solid line is a fit result of a Gaussian peak plus a background function. The residual background is shown as a dotted 
fine. Panels (a) and (b) are from a low pr bin in MB events; panels (c) and (d) are from a high pr bin in HT events. Panels 
(b) and (d) are distributions after mixed-event background subtraction from panels (a) and (c). 



An EMC photon is reconstructed from a single tower. 
The photon energy is determined by the tower energy. In 
MB and central events, towers with energy greater than 
500 MeV and at least 250 MeV higher than any of their 
eight surrounding towers are required. The BSMD hit in- 
formation is not used due to its expected inefficiency for 
low energy photons. The photon position is assumed to 
be at the center of the tower. Charged particle contami- 
nation is greatly reduced by projecting TPC tracks into 
the BEMC and vetoing the first two towers intersected 
by the track. 

In HT triggered events, BSMD hits are used to sepa- 
rate the two close decay photons from a single 7r° decay. 
The photon positions are determined from the BSMD 
hits. If multiple photons are found in the same tower, the 
tower energy is split according to the individual BSMD 
hit energies. For a photon with energy below the HT 
threshold, we require that no TPC track is projected into 
an area of ± 0.05 in Ar/ and ± 0.05 in A(j> around the 
photon candidate. For a photon above the HT threshold, 
we require that the sum of momenta of all charged par- 
ticle tracks projected to the surrounding Ar] — A</> area 
should be less than 1 GeV/c. 

For TPC photons, we select e + /e~ candidates via ion- 
ization energy loss dE/dx in the TPC. A number of ge- 
ometrical cuts are applied to each e + /e~ pair to have a 
topological signature of a photon conversion. These cuts 
require that the two tracks originate from a common sec- 
ondary vertex within or before entering the TPC with a 
small opening angle and a small invariant mass, and that 
the reconstructed photon candidate originates from the 
primary vertex. The photon momentum is taken as the 
sum of two daughter track momenta at the conversion 
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FIG. 2: The di-photon invariant mass distributions using the 
EMC-EMC method in 0-20% Au+Au collisions at ^/snn = 
200 GeV. The solid line is a fit result of a Gaussian peak plus 
a background function. The residual background is shown as 
a dotted line. Panels (b) is the distribution after mixed-event 
background subtraction from panel (a). 



point. This technique has been used in Au+Au collisions 
at y / SNN = 130 GeV, and more details can be found in 
Rcf. [10]. 
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7r Reconstruction 



In MB and central events an EMC photon is paired 
with a TPC photon (EMC-TPC), and in HT events pairs 
of two EMC photons (EMC-EMC) are also used to re- 
construct 7r°'s. Due to the relatively large tower size of 
the STAR BEMC, energy determinations for low px pho- 
tons can be contaminated due to the presence of other 
nearby photons, neutral particles, and residual charged 
particles in high multiplicity Au+Au events. The limited 
energy resolution of BEMC towers at low energy (nom- 
inal resolution of BEMC towers has been estimated to 
be 16%/ \f~E © 1.5%) [8] also hinders the accurate mea- 
surement of photon energy. As a result, it is difficult to 
obtain a clear 7r° signal at low px by exclusively pairing 
EMC photons. On the other hand, the relatively tight ge- 
ometrical cuts for TPC photon reconstruction select very 
clean conversion photon samples. They significantly re- 
duce the combinatoric background and improve the ir° 
mass resolution. The EMC-TPC method yields a clear 
7T° signal from 1 GeV/c to intermediate px (~5 GeV/c) 
in central Au+Au collisions. At higher px above the HT 
threshold, the EMC photons are less affected by back- 
grounds and the EMC-EMC method produces clear 7T° 
signals. Due to its greater efficiency for high px photons, 
the EMC-EMC method is able to extend the measure- 
ment to higher px ■ 

The mixed-event technique is used to reproduce a com- 
binatoric background. For the mixed-event distribution, 
photons from the event being analyzed are paired with 
photons from events in an event pool, in which events are 
required to have similar multiplicity and primary vertex 
position as the one being analyzed. The di-photon in- 
variant mass distribution after mixed-event background 
subtraction is fit to extract the raw ir° yield. 

Figures 1 and 2 show examples of the di-photon in- 
variant mass distributions before and after mixed-event 
background subtraction for different ir° reconstruction 
methods. In Fig. 1 we show the invariant mass distribu- 
tions in the px regions 1.2-1.5 and 7.0-9.0 GeV/c from 
the EMC-TPC method, and Fig. 2 shows the invariant 
mass distributions in the px region 9.0-12.0 GeV/c from 
the EMC-EMC method. All invariant mass distributions 
are for the 0-20% collision centrality bin. These distri- 
butions are fit using a Gaussian plus a polynomial back- 
ground function. At highpr, the background is small and 
can be easily subtracted by fitting a linear dependence 
on Mi nv . At low px, the signal-to-background ratio is 
rather small. After mixed-event background subtraction, 
a larger residual background is observed and a 3rd order 
polynomial function is used to fit the background shape. 
Here the normalization factor between same-event and 
mixed-event is adjusted for eachp-r bin, so that the resid- 
ual background has a shape roughly linearly increasing 
with mass, and can be described by a polynomial fit. The 
residual background may come from correlated photons 
that are not combinatoric and cannot be reproduced by 
the mixed-event technique. Such correlations may arise 
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FIG. 3: The 7r° invariant mass peak positions (a) and peak 
widths (b) as a function of pr in 0-80% Au+Au collisions at 
^/sjvjv = 200 GeV. Corresponding results from 7r° embedded 
simulation are shown for comparison. 



from contaminations to EMC photons or from resonance 
decays to multiple photons in the final state. The amount 
of residual background is strongly centrality and px de- 
pendent, more pronounced in the most central events and 
at lower px- Figure la shows the situation where both 
combinatoric and residual backgrounds are most severe; 
nevertheless the ir° signal can still be observed above the 
residual background (Fig. lb). 

Several systematic checks have been performed. 
Firstly, a track rotation technique was used to gener- 
ate combinatoric background for comparison. It rotates 
the EMC photons by 180° in the azimuthal plane, and 
mixes them with the photons reconstructed in the TPC. 
Secondly, the normalization factor was adjusted and the 
invariant mass distribution was re-fit to extract ir yield. 
Although these two procedures may significantly change 
the shape of residual background, yields extracted using 
the same function are consistent with each other. We 
have also changed the order of polynomial used for back- 
ground fitting, as well as the fit range, and have included 
the variance in the overall systematic errors. 

Figure 3 shows the extracted 7r° peak positions and 
widths as a function of px using different 7r° reconstruc- 
tion methods. Results from real data are compared to 
Monte Carlo (MC) simulation of w° embedded in real 
data. The n° peak position shows some px dependence 
at low px for the EMC-TPC method and at higher px 
for the EMC-EMC method. At low p T the drop is un- 
derstood as the effect of energy loss of e + /e~ tracks due 
to bremsstrahlung. At high px the rise of peak position 
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FIG. 4: Left: Photon conversion point radius distributions from real data and MC simulation in Au+Au collisions at ^snn 
— 200 GeV. The two distributions are normalized in the radius region of the TPC gas at 60 < r < 100 cm. Right: Conversion 
probability correction factor for n° as a function of pr- 



os a function of pr is due to the saturation of dynamic 
range for energy measurement from the BSMD in this 
data set. The BSMD read-out saturated when the de- 
posited energy exceeds about 6 GeV, leading to more 
evenly distributed energies of two spatially close photons 
when they hit the same BEMC tower, and therefore pro- 
duces a larger invariant mass. This effect is more pro- 
nounced at higher pr- The saturation scale was lower 
than anticipated due to electronics signal termination is- 
sues in 2004. The effect has been included in the simu- 
lation. The trend of pt dependence is well reproduced 
by the simulation but the simulation underestimates the 
mass peak position by 4-8%. The use of TPC photons 
significantly improves the ir° peak resolution. For the 
same HT data sample the EMC-TPC method yields peak 
widths narrower than those from the EMC-EMC method, 
which is consistent with the MC simulation. Comparing 
the MB and HT data samples, the requirement of BSMD 
hits improves the spatial resolution of EMC photons, and 
thus measures a narrower ir° peak width. 
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FIG. 5: The overall detection efficiency of n a from embedding 
study in 0-80% Au+Au collisions at ^sj^ = 200 GeV. 



D. 7T° Detection Efficiency 

The raw yield of 7r° is corrected for an overall central- 
ity dependent detection efficiency calculated from a full 
MC simulation, embedding 7r°'s into real events. The 
efficiency is calculated as the ratio of the reconstructed 
ir° pt spectrum over the input spectrum, using the same 
cuts as the real data analysis. The input ir° spectrum 
for the embedding analysis is weighted so that it repro- 
duces the previously measured charged [12] or neutral 
pion spectrum [13]. The calculated efficiency takes into 
account the losses due to acceptance, photon conversion 
probability, tracking inefficiency, track energy loss, and 
track quality cuts. 

The conversion probability is crucial for the ir° effi- 



ciency calculation when TPC photons are used. A com- 
parison of distributions for the photon conversion radii 
between data and MC simulation is shown in the left 
panel of Fig. 4. The two distributions are normalized to 
the data in the inner radius region of the TPC gas where 
geometry is simple and well understood. The comparison 
indicates that the photon conversion probability can be 
well reproduced in the regions of TPC gas and the inner 
field cage, but is underestimated in the regions of SVT 
and SSD where structures are complex. Similar observa- 
tion has been reported in earlier publications [10]. The 
results for TPC photons are corrected for this effect, us- 
ing the conversion rates in the TPC g clS 0jS el reference. A 
correction factor F geo = {n det /n gas ) da ta/ '(n det /n gas ) MC 
is calculated as a function of the conversion photon pt, 
where n de t and n gas are numbers of conversion points 
in the whole detector and in the TPC gas only. In the 
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embedding analysis a reconstructed TPC photon associ- 
ated with a MC photon is weighted by the factor F geo 
corresponding to its px, which folds the correction in the 
efficiency calculation. The final correction factor for tt° 
as a function of px is shown in the right panel of Fig. 4. 

Figure 5 shows the overall detection efficiency as a 
function of px in Au+Au collisions at ^/snn = 200 
GeV. The use of TPC photons is statistically challenging 
due to the relatively low conversion probability. Using 
EMC photons enhances the efficiency significantly and is 
preferable in studying the ir° spectrum at high px- The 
efficiencies shown here have taken the conversion proba- 
bility correction into account. 



E. Systematic Errors 



Major sources of systematic errors for the 7r° mea- 
surement in Au+Au collisions at ^Jsnn = 200 GeV are 
listed in Table I. Systematic errors are calculated for each 
px bin, and systematic errors from different sources are 
added in quadrature. The systematic errors are esti- 
mated by using several methods. First, we have varied 
photon reconstruction cuts and compared the fully cor- 
rected spectra. By changing the geometrical cuts applied 
to the TPC photon reconstruction, a systematic error of 
10-20% in the final spectra is obtained. For the EMC 
photons in MB events, different energy cuts are used to 
select EMC photon samples with different levels of pu- 
rity. In HT events, instead of a single tower, a cluster 
with up to 2 x 2 towers is used to reconstruct an EMC 
photon. These various photon reconstruction methods 
give a systematic error of 10-20% to the final tt° spectra. 
Next, We have also included the uncertainties in the ab- 
solute energy scale of the BEMC, which could affect the 
overall shape of the ir° spectra. An offset in the BEMC 
energy scale would contribute to the small deviations in 
the 7T° mass peaks between real data and simulations. We 
have included a uncertainty of ±5% [14] in the BEMC 
energy scale in our Monte Carlo simulations. We esti- 
mated systematic Errors of 20-35% throughout the px 
range. Third, we have varied the fitting procedure used 
to extract the raw ir° yield. The raw yield of 7r° depends 
on the background fitting function, fit range, and the 
normalization of same- and mixed-event invariant mass 
distributions. Results using different fitting parameters 
indicate a systematic error of 10-15%. We have also cross 
checked the uncertainty due to the conversion probability 
correction by applying the correction factor as a function 
of conversion point position. The result agrees with the 
original within 10%. As a result, a 10% systematic er- 
ror is assigned for the conversion probability correction 
factor. 



Q. 
■D 



Q. 

N 
CM 



10 

1 

10" 2 
10" 4 
10' 5 
10" 8 
lO" 10 
10" 12 



-I 


I 




I 


I I - 




*> 




Au+Au at\[s^ 


= 200 GeV 




AO 




o A STAR MB EMC-TPCji 




A O 


• a STAR HT EMC-TPCji ~ 


- 


oV 


STAR HT EMC-EMCtc ~ 




A°D A O 


0-1 2% x 10 1 






A O 


















^ A f 








O 










o i 


* ^ ir 


MB (0-80%) - 






t 


t * 












0-20% x 10" 2 1 








# 


20-40% x 10" 3 J 










40-80% x 10 -4 - 





2 


4 


6 8 10 12 14 16 18 



Transverse Momentum p T (GeV/c) 



FIG. 6: Invariant yield of STAR tt° as a function of pr at 
mid-rapidity for different collision centrality bins in Au+Au 
collisions at ^/snn = 200 GeV. Spectra for different colli- 
sion centralities are scaled for clarity. Statistical errors are 
shown as vertical lines and point-to-point systematic errors 
are shown as bars. 



III. RESULTS 

The 7T° invariant yield per collision at mid-rapidity 
(0 < y < 0.8) as a function of px in Au+Au collisions 
at y'SNN = 200 GeV is shown in Fig. 6. Statistical and 
systematic errors are shown as vertical lines and bars, 
respectively. The horizontal size of the vertical bars also 
indicates the range of the px bin. In addition to the 
overall MB 0-80% result, the data sample is also divided 
into three collision centrality bins 0-20%, 20-40%, and 
40-80% based on measured charged particle multiplicity 
at mid-rapidity from the TPC [15], with 0-20% the most 
central collisions. The ir° spectra are measured over an 
extended px range from 1 to 12 GeV/c. Results from 
different ir° reconstruction algorithms and different data 
samples were compared in overlapping px ranges, and 
were found to be in good agreement. Therefore, in the 
following figures only a combined data point using sta- 
tistical weighted average of data from various algorithms 
will be shown in the overlapping px bins. 

Our 7T° spectra are compared to the previously pub- 
lished tt^~ and 7T° results. The ratios of our measured 7r° 
spectra to the STAR tt* [12] and the PHENIX tt [13] 
in Au+Au collisions at ^/snn = 200 GeV are shown in 
Fig. 7. Parametrized results of the tt° and 7r ± spectra 
from power law functions are used in order to match the 
px binning of our ir° data. The error bars are propa- 
gated using the averaged error of two neighboring data 
points. The spectrum ratio is slightly larger in periph- 
eral collisions than in central and mid-central collisions. 
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TABLE I: Summary of main sources of systematic uncertainties on the 7r° yields in Au+Au collisions at y^ijvjv = 200 GeV. 
Systematic uncertainties from varying photon cuts and yield extraction techniques are pr uncorrelated, and systematic uncer- 
tainties from the BEMC energy scale and conversion probability correction are pr correlated. 
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FIG. 7: The ratios of STAR n° spectra over tt* from STAR [12] (solid symbol) and it from PHENIX [13] (open symbol) as a 
function of pr for different collision centrality bins in Au+Au collisions at y / ijvjv = 200 GeV. Statistical errors are shown as 
vertical lines and point-to-point systematic errors are shown as bands. 



With the best statistics in MB 0-80% collision central- 
ity, the STAR tt° yields are consistent with the PHENIX 
7T° yields, and about 15% smaller than the STAR 
yields over the pt range. Considering that the system- 
atic uncertainties in the STAR ir° and ir^ analyses are 
mostly independent, the two yields are consistent within 
systematic uncertainties. 

The nuclear modification factors can be calculated us- 
ing peripheral collisions as a reference (Rcp) or using 
nucleon-nucleon collisions as a reference (Raa)- 



1 central 



Rcp{pt) 



and 



Raa(pt) 



[d 2 N/ PT dyd PT /(N hin ) 
[d 2 N/p T dydp T /(N hin )]P e "P heral ' 

d 2 N AA /dydp T /{N hin ) 



where (iVbin) is the average number of binary nucleon- 
nucleon collisions per nucleus-nucleus collision. The a. 



d 2 a pp /dydp T /ui 



VP 

The 



is taken to be 42 mb for ^sii = 200 GeV [16]. 
measurements of suppression for high px charged hadrons 
from STAR [1] and neutral pions from PHENIX [2] (R CP 
and Raa < 1) in most central Au+Au collisions at RHIC 
provided the first experimental evidence that partons suf- 
fer energy loss in the dense matter created in these colli- 
sions. 

Figure 8 shows our measurements of the nuclear mod- 
ification factor Rep for ir° as a function of pr in Au+Au 
collisions at ^s^v = 200 GeV for the 0-20% and 20-40% 
over 40-80% collision centrality bins. When calculating 
Rcp some systematic uncertainties cancel out, such as 
the BEMC energy scale and conversion probability cor- 
rection. Compared to the 40-80% peripheral Au+Au col- 
lisions, the more central collisions show a suppression of 
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FIG. 8: The nuclear modification factor Rcp as a function 
of p T of STAR ty° compared to STAR vr ± [12] in Au+Au 
collisions at y / ijvjv = 200 GeV. Statistical errors are shown as 
vertical lines and point-to-point systematic errors are shown 
as solid lines. The shaded band on the right demonstrates 
the uncertainty of iVbin- The dashed curves are jet quenching 
theoretical calculations [18]. 



the 7T° yield indicated by Rcp < 1 and the suppression 
is even stronger for the most central collisions. At high 
Pt > 4 GeV/c the ir° Rcp is independent of pt within 
uncertainties. Our measured 7T° Rcp values show the 
same magnitude of suppression as the STAR tt ± data 
[12], which are shown as open circles. 

Figure 9 shows our measurements of the nuclear mod- 
ification factor i?AA for ir° as a function of pt in Au+Au 
collisions at v / s A 7 A 7 = 200 GeV for 0-20%, 20-40%, and 
40-80% collision centrality bins, where a parameterized 
description of the tt^ spectrum in p+p collisions at ^/snn 
= 200 GeV from Ref. [17] is used to calculate R AA . The 
error bars of the p + p result are propagated using the 
averaged error of two neighboring data points. The ir° 
Raa shows a similar centrality dependence as the Rcp- 
In the most central Au+Au collisions the ir° yield is sup- 
pressed by a factor of about 5 relative to the expectation 
from scaled nucleon-nucleon collisions. For all the colli- 
sion centrality bins, our measured Raa values for ir° at 
high pt agree with previously published results from the 
PHENIX collaboration [13] within systematic uncertain- 
ties. 

The nuclear modification factors for inclusive light 
quark mesons at high pt in central heavy-ion collisions 
have been investigated with several model calculations. 
The nuclear modification factors Rqp and Raa have been 
calculated in terms of parameters such as the initial gluon 
density [18, 19] and the medium transport coefficient q 
[20], which characterize properties of the dense matter 
created. Recent theoretical calculations suggest that the 
collisional energy loss may also play an important role 



2 P 

I * STAR 7i° 0-2"% 

5 - O PHENIX 7i° 

I Vitev R.. 0-10% 



< 
< 
DC 




5 10 15 

Transverse Momentum p t (GeV/c) 



FIG. 9: The nuclear modification factor Raa as a function 
of p T of STAR 7T° compared to PHENIX tv° [13] in Au+Au 
collisions at y^ijvjv = 200 GeV. Statistical errors are shown as 
vertical lines and point-to-point systematic errors are shown 
as solid lines. The shaded band on the right demonstrates 
the uncertainties of iVbin and the normalization uncertainty 
in p+p collisions of 14% [3] . The dashed curves are theoretical 
calculations in 0-10% Au+Au collisions [19]. 



in explaining the large suppression of non-photonic elec- 
trons from heavy quark decays [21]. In Fig. 8, we show 
an example of a theoretical calculation of Rqp with ini- 
tial gluon density dN 9 /dy = 1150 in 0-10% Au+Au and 
between 100 and 150 in 40-80% Au+Au collisions [18]. 
In Fig. 9, theoretical calculations with dN 9 / dy— 800 
to 1150 for 0-10% Au+Au collision centrality [19] are 
shown as dashed curves in comparison to measurements 
from STAR and PHENIX from 0-20% collision centrality. 
Experimental measurements and theoretical predictions 
agree reasonably well, indicating that the yield suppres- 
sion of light quark mesons may be accounted for by the 
parton energy loss mostly through gluon radiation. 



IV. SUMMARY 

We have presented the first STAR results for 7r° pro- 
duction in Au+Au collisions at ^/snn = 200 GeV. The ir° 
spectra are measured over the range of 1< pt <12 GeV/c 
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using the combination of conversion photons from TPC 
reconstruction and photons from BEMC energy measure- 
ment. Despite the relatively large tower size, the STAR 
BEMC alone can be used to reconstruct 7r°'s for px > 4 
GeV/c. The use of conversion photons significantly en- 
hances detection capability for 7r°'s at low and interme- 
diate pt and extends our tt° measurement to a lower px 
range. 

Our measurements of 7r° spectra are consistent with 
the STAR charged tt* and PHENIX tt° results within 
statistical and systematic errors. The nuclear modifica- 
tion factors Rcp and Raa of the STAR 7r° data confirm 
the previously published ir results and can be described 
by theoretical calculations based on parton energy loss 
through gluon radiation in the dense medium created at 
RHIC. In the most central Au+Au collisions the inclu- 
sive 7T° yield shows a factor of about 5 suppression rel- 
ative to the expectation from scaled p + p collisions for 
Pt > 5 GeV/c. Our measurements confirm the mag- 
nitude of light hadron suppression observed in central 
Au+Au collisions and provide further support for the 
physical picture of jet quenching in the dense matter cre- 



ated in nucleus- nucleus collisions at RHIC. 
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